Quantum information processing using localized ensembles of nuclear spins 



O 

o 



(N 



c3 



-a 
o 

(N 
> 
O 

o 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



J. M. Taylor 1 , G. Giedke 2 , H. Christ 3 , B. Paredes 3 , J. I. Cirac 3 , P. Zoller 4 , M. D. Lukin 1 , and A. Imamoglu 2 

1 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 
2 Institut fur Quantenelektronik, ETH Zurich, Wolfgang- Pauli-Strafle 16, CH-8093 Zurich, Switzerland 

3 Max- Planck- Institut fur Quantenoptik, Hans-Kopfermann-Strafie 1, Garching, D-85748, Germany 
4 Institut fur Theoretische Physik, 



Wolfgang-Pauli-Strafie 16, 
;s-Kopfermann-Strafle 1, G 
Universitat Innsbruck, Technikerstrafie 24, A-6020 Innsbruck, Austria 




We describe a technique for quantum information processing based on localized ensembles of 
nuclear spins. A qubit is identified as the presence or absence of a collective excitation of a mesoscopic 
ensemble of nuclear spins surrounding a single quantum dot. All single and two-qubit operations 
can be effected using hyperfine interactions and single-electron spin rotations, hence the proposed 
scheme avoids gate errors arising from entanglement between spin and orbital degrees of freedom. 
Ultra-long coherence times of nuclear spins suggest that this scheme could be particularly well suited 
for applications where long lived memory is essential. 

PACS numbers: 03.67.Lx, 71.70Jp, 73.21. La, 76.70.-r 



Nuclear spin degrees of freedom have attracted consid- 
erable attention as potential carriers of quantum infor- 
mation due to their exceptionally long coherence times. 
Early bulk NMR work [lj has substantially enriched our 
understanding of the key features of quantum computa- 
tion [2|, |3( • The fundamental difficulties in scaling bulk 
NMR to a large number of qubits motivated efforts to 
use single, individually addressable nuclear spins in semi- 
conductors as qubits , where computation is primarily 
mediated by the hyperfine interaction between electron 
and nuclear spin. While possibly scalable, such a scheme 
is limited by the fact that the electron wave-function is 
spread over many lattice sites, reducing the strength of 
the hyperfine interaction. In addition, two-qubit oper- 
ations in Ref. rely upon exchange coupling, making 
them susceptible to fast orbital decoherence mechanisms. 

Recently, a method for robust storage of quantum in- 
formation in localized ensembles of nuclear spins was sug- 
gested 0, , where it was shown that the collective hy- 
perfine coupling between nuclear and electron spin de- 
grees of freedom provides a controllable mechanism for 
coherent storage and manipulation of quantum states. 
These nuclear spin ensembles correspond, for example, to 
the lattice nuclei in a quantum dot. As a quantum mem- 
ory, such nuclear ensembles are robust with respect to 
variations in dot characteristics, rely upon proven fabri- 
cation techniques, and provide high fidelity storage with- 
out requiring a high-degree of nuclear spin polarization. 

In this Letter, we describe a technique to efficiently 
process quantum information stored in localized nuclear 
spin ensembles. Specifically, these ensembles enable a ro- 
bust, scalable implementation of quantum computation 
protocols, unencumbered by the difficulties faced by sin- 
gle spin impurity or bulk NMR approaches. The fun- 
damental interaction that allows for spin manipulation 
in our scheme is hyperfine coupling: as a result, the or- 
bital and spin degrees of freedom remain unentangled 
throughout the two-qubit gate operations, mitigating the 
effects of orbital decoherence on gate fidelity. While col- 



FIG. 1: Quantum dot based approaches for nuclear ensemble- 
based computation, (a) Electrically defined (lateral) quantum 
dots connected by ballistic transport e.g. through a series of 
quantum dots. Nearby rf-SETs (not shown) would provide 
measurement, (b) Nanowhiskers Q with several optically ac- 
tive dots. Tasks between A and A' are achieved via shuttling 
of an electron from the reservoir dot via a potential difference, 
Va- Long distances tasks are performed through photon- 
based entanglement generation between A' and B'. |19j 



lective enhancement of hyperfine interaction allows for 
fast quantum gates, the ultra-long nuclear spin coher- 
ence times render the scheme particularly attractive for 
memory intensive quantum information processing tasks. 

The scheme proposed here can be realized using either 
electrically || or optically manipulated quantum dots 
that are defined as in Fig. ^ The nuclear ensemble in 
each dot is prepared using polarized electron spins |rjflfl|. 
We illustrate our technique by considering the fully po- 
larized case (when the nuclear state is |0) = |— I, . . . , —I) 
for N spin-/ nuclei), demonstrate how to perform single 
qubit and two-qubit gates, and consider sources of error. 

The Hamiltonian describing a single electron interact- 
ing with the nuclear spins in a quantum dot is |I1| 

H = hB ls S z + hB li k) 4 k) + ha S -A. (1) 



The first and second terms describe the coupling of the 



electron and the nuclear spins to an external magnetic 

(k) 

field, with effective gyromagnetic ratios 7s, 7} . The 
last term is the electron coupling to collective nuclear 
degrees of freedom, defined by A±,z = J2k a kl±\- The 
a/c's are proportional to the probability density of the 
electron at the location of the corresponding nuclei and 
are normalized such that J2k a k — N; a — A/N is the 
average hyperfine interaction per nucleus. Both A and N 
depend upon the specific material and dot construction; 
typical numbers are A ~ 10-100 ns _1 and N = 10 4 -10 6 
nuclei. Excitations out of the fully polarized state form 
an orthonormal set of collective nuclear spin states |m) oc 
(A + ) m |0), where m is the number of excitations. 

When an electron is confined in the dot, evolution over 
times much shorter than a -1 is restricted to subspaces 
spanned by {|m — 1) ||) , \m) |J,)}. Using Pauli matrices, 
the Hamiltonian for each subspace with a fixed excitation 
number m is H^ m ^ = h5 m a z m ^ + ftQ, m <j£"^ , with 

ti m = aj2\j (m - 1| \m - 1); (2a) 

6 m = a/i{(m\A z \m) + {m-l\A z \m- 1)) (2b) 
+ [j s ~ (H K z \m) - (to - 1| K z \m - l))]B/2, 

where K z — ^2"/j k ^I z k ^ is a sum over nuclear spin op- 
erators, weighted by individual nuclear spin gyromag- 
netic ratios. When Overhauser shift and Zccman energy 
sum to zero (\S m \ <C fl m ) and coherent flip-flop (Rabi) 
oscillations occur at rate O m . The energy level struc- 
ture of the coupled electron-nuclear system is shown in 
Fig.|21along with the coupling strengths for m < 2: since 
&m+i = Vm^-m where i] m = y/m + 1[1 - 0(m/N)], it 
is easy to note the analogy with the celebrated Jaynes- 
Cummings (JC) model of quantum optics 01 • We use 
the nonlinearity of such a JC-type two-level system cou- 
pled to a nearly-bosonic mode to effect elementary quan- 
tum gates. 

Quantum information stored in the m = 0, 1 manifold 
can be mapped reliably from nuclear states to electron 
spin and back via a generalized rotation : 

^(tt/2, 0)(a |0) + 11)) ||) = |0) (a ||) + |T», (3) 

where 

i£»(0, 4>) = e i*5* e -<flH/(«ii) e -i*S. . (4) 

The transfer of quantum information from the nuclear 
ensemble to electron spin allows for fast single qubit 
operations to be performed: after a \\) electron is in- 
jected into the dot, the quantum information is trans- 
ferred to the electron spin, and then the operation is 
performed on the electron. Finally, the quantum infor- 
mation is mapped back to the nuclear ensemble. A z- 
axis rotation can be accomplished by waiting in a static 
magnetic field or through a laser-induced spin-dependent 
AC Stark shift M x-axis rotations can be done via 
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FIG. 2: Level structure of the combined electron spin and nu- 
clear dark state plus excited states. For N — * 00, the coupling 
between the two excitation manifold is stronger by -q — > v2. 

ESR[i[3 (with ft ES R> 1 

ns 1 ) or optical Raman spin 
flips through a virtual trion state (with f2 pt.ESR ^ 10 
ns_1 ) 01- Measurement of the ensemble nuclear spin 
state can be implemented by mapping the quantum in- 
formation to electron spin, and carrying out an electron 
spin measurement either by state selective ionization fol- 
lowed by charge measurement with a rf-SET |14| or by 
detect ing fluorescence in a spin-dependent cycling tran- 
sition [miii]. 

To perform a two-qubit gate between quantum dots A 
and A', a single electron can be used as to transfer quan- 
tum information between the dots: the state of nuclear 
spin qubit A is mapped onto the electron, which is then 
moved to A', where a controlled-phase (CP) gate between 
the nuclear and the electronic qubit is applied using the 
nonlinearity of the interaction. Following Ref. [ljj, a 
two-qubit CP-gate (up to single qubit gates) is given by 

^(tt/4, mtl^h, -7r/2)iC(-V4,0). (5) 
In the computational basis, this corresponds to 

/l \ 

e i7r /" . . 

e-™l ri ■ ^ ' 

\0 -\j 

After this operation, the quantum information carried by 
the electron spin is mapped back onto the nuclear spin 
ensemble A. Applying a z-rotation of —tt/t] to qubit A 
and one of n/r] to qubit A' yields a CP gate. 

For distant dots, where electron spin transport be- 
tween ensembles is difficult or impossible, it may be 
easier to generate entanglement off-line, apply local pu- 
rification protocols and then use it to effect non-local 
gates, following 18]. Entanglement between distant elec- 
tron spins can be generated optically, using spin-flip op- 
tical Raman transitions [lflj . or electronically, such as 
through adiabatic splitting of a spin singlet in a double 
dot |20j. The electron- nuclear state mapping procedure 
can then ensure that the distant nuclear spin ensembles 
are entangled. Starting with the entangled nuclear state 
(|01) + 1 10>) (between dots A' and B', cf. Fig. [TJ, we 
can implement a CP gate on dots A and B (determinis- 
tically) by performing local unitaries and measurements 
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of AA' and BB' as follows: (i) perform a CNOT A ^A> 
then measure A', (ii) Perform a CP gate 1 10} — > — 1 10) 
at BB' followed by a Hadamard gate and measurement 
at B'. (iii) Local phase flips |0) -> - |0) at A (B) if the 
measurement outcomes at B'(A') were "1" complete the 
CP gate between A and B [2l|. 

We next analyze various sources of errors due to finite 
polarization, inhomogeneity, nuclear spin dynamics, and 
electron spin decoherence. To understand the role of fi- 
nite polarization, the specific cooling procedure must be 
examined. When the nuclear ensemble starts as a ther- 
mal mixture, cooling to dark states can be achieved by 
coupling polarized electron spins to the nuclear ensem- 
ble [lfj. Regardless of the details of cooling, the final 
density matrix will be a statistical mixture of dark states 
\T>(n,f3)), where A_ \T>(n,/3)} = 0. n is the number of 
spin excitations (n — is the fully polarized state) and (3 
is the permutation group quantum number |2^ . It was 
recently shown 6] that dark states have the same symme- 
try properties as the fully polarized state, and a manifold 
of excited states can be defined from a given dark state, 
\m(n,(3)) oc (A + ) m \T>(n,f3)). Hence the above consider- 
ations for perfectly polarized nuclei map directly to the 
case when the nuclear ensembles start in any given dark 
state, not just the fully polarized state. 

In practice, the cooled nuclear ensemble density ma- 
trix is a mixture of different dark states, i.e. p = 
J2 n a Pn,p\'D(n,f3))('D(n,/3)\. As each dark state has a 
different d m (n, /3),£l m (n, P), interaction times and ap- 
plied magnetic fields can only produce R^($,cf)) with 
the desired angle 9 for some fraction of the given mix- 
ture. The inhomogeneous mixture effects can be charac- 
terized by examining the subgroups of dark states with 
different detunings, which lead to errors in Rabi oscilla- 
tions p ~ (cra/fii) 2 , where as is the standard deviation 
of possible 5% values over the distribution p n .p- In the 
homogeneous case with spin- 1/2 nuclei, 

a 5 ~ a^N(l - P 2 ) . (7) 

Even at high (P ~ 0.95) polarizations, the effect of dif- 
ferent detunings can be substantial (p ~ 0.03). This 
provides the strongest limit to realization of 23J . 

The inhomogeneous nature of the hyperfine coupling 
leads to further errors. In this case, the logical states 
of the system (the m-excitation manifolds) are no longer 
eigenstates of J z = A z +S z and of A 2 . As a consequence, 
there is a nonzero probability that the system moves out 
of the computational subspace during the gate opera- 
tion. We estimate these leakage errors using the tech- 
niques developed in Refs. 0, @ , and find that the result- 
ing gate error Pmhom decreases with increasing number 
of nuclei: for N ~ 10 5 at high polarizations (P > 0.95) 
Pinhom ^ 10~ 3 [lil . We note that a similar error emerges 
due to the differences in Zeeman energy associated to dif- 
ferent nuclear species. For materials like GaAs, with gy- 
romagnetic ratios varying greatly from species to species, 



this limits the effectiveness of gate operations at high 
magnetic field, resulting in the errors in the range of 
10~ 3 as indicated in Fig. 3. Optical manipulation (e.g. 
tuning the system into resonance via spin-dependent AC 
optical Stark shifts) may mitigate this difficulty. Finally, 
in between gate operations the errors associated with in- 
homogeneous evolution may be eliminated by refocusing 
sequences of NMR pulses. 

The nuclear spin diffusion due to dipole-dipole interac- 
tions with rate jdd ~ 60ms -1 leads to a decay of the co- 
herences which form the dark states [2j| . Active correc- 
tions can be performed with NMR pulse sequences that 
average the dipole-dipole Hamiltonian, such as WaHuHa, 
improving the error rate to T 2 hh j^ D |26j,|2jj. After cor- 
rection dark state coherences could have lifetimes on the 
order of 0.1-ls for moderate cycle times r w hh- Quadrupo- 
lar terms due to inhomogeneous strain can cause ad- 
ditional differential phase evolution of different nuclear 
spins. However, this type of phase inhomogeneity leads 
to an error second order in the interaction which is neg- 
ligible (~ 10~ 7 per cycle of computation). 

The errors of manipulation of single electron spins in 
quantum dots via microwave or optical fields have been 
considered in detail elsewhere 0, 0, El El , and only the 
relevant results are quoted here. Electron spin decoher- 
ence will most likely be limited by different Overhauser 
shifts corresponding to different detunings, with error go- 
ing as p ~ (ct^/^esr) 2 ; fast ESR (~ 6ns~ 4 ) will mitigate 
this effect; with optical fields, even faster effective Rabi- 
oscillations are predicted (Oopt.ESR ~ 50ns -1 ) with er- 
rors then limited by spontaneous emission to ~ 10~ 3 . As 
for measurement, the fundamental limit will be set by 
relaxation of the electron spin, which has a time scale 
r _1 ~ 0.1-10ms, based on recent measurements of the 
spin relaxation time 0, . 

Moving individual electron spins over short distances 
incurs phase errors due to uncertainty in the local nuclear 
field of the transport channel. The error induced by ran- 
domly oriented nuclei during transit can be estimated by 
the time-averaged hyperfine field the electron wavepacket 
encounters over the whole process, A/y nLhl, where n is 
the density of lattice nuclei, L is the length of the chan- 
nel, and hi is the transverse area of the channel. For a 
50nm wide, lOnm high channel, and a quantum dot sep- 
aration of 1/im, the expected dephasing probability due 
to thermal nuclei is ~ 5 x 10 -5 for a 3ns transfer time. 

We consider three materials with demonstrated elec- 
tronic and optical quantum dots (GaAs, N = 10 5 ; InAs, 
N = 10 4 - 6 ; CdSe, N = 10 4 ). The expected error of a 
two-qubit gate operation, defined as Tr[J7 pcr f 0Ct J7^ ctual ]/4, 
where f7 pe rfect is given by Eqn.El is plotted in Fig. 03 As 
the error is dominated by detuning error, it is material 
independent, and is a few percent at 95% polarizations. 

To combat these errors, a series of measurements may 
be made to determine the effective detuning of the sys- 
tem better than the thermal mixture limit as- This can 
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FIG. 3: For a cooled dark state density matrix (c.f. 
expected error for a controlled-phase operation between elec- 
tron spin and nuclear spin as a function of polarization (CdSe: 
solid, InAs: dotted, GaAs: dashed; in the main figure, CdSe 
and InAs are indistinguishable). Inset shows effects of nar- 
rowing by a factor of 10 for different materials. 

be done e.g. by a measurement made with the single elec- 
tron in the quantum dot, in direct analo gy t o single ion 
Ramsey interferometry in atomic clocks |30| : the distri- 
bution is then narrowed, with a$ ~ as I \fn, where n is the 
number of measurements. An improvement in the uncer- 
tainty of the mixture's detuning of a factor of 10 yields 
high fidelity operation, as shown in the inset to Fig.|3 31] . 
Smaller dots produce better results at higher polariza- 
tions due to their greater coupling strengths. The limits 
for GaAs are due to the large species inhomogeneity and 
the incommensurate requirement of sufficient magnetic 
field to perform effective coherent averaging of dipole- 
dipole interactions. Low species inhomogeneity allows for 
higher magnetic field, faster WaHuHa-type correction se- 
quences, and fewer errors in the transfer operation. This 
would be the case for quantum dots defined in nanotubes 
with isotopically enhanced 13 C or in CdSe quantum dots. 
Materials with low spin-orbit interaction will also reduce 
electron spin dephasing. 

In conclusion, we have detailed a scheme for quantum 
information processing using dynamically defined qubits 
composed of collective excitations of nuclear spins in a 
quantum dot. For small dots with near- homogeneous 
Zeeman splittings (e.g. CdSe, InAs) large but finite po- 
larization (95%) is already sufficient to reach an error 
rate in two-qubit operations of order a few percent. 
Given the suppression of dipolar diffusion in lithograph- 
ically isolated structures (e.g. vertical quantum dots, 
nanowhiskers, self-assembled quantum dots) these po- 
larizations may be within reach. Moderate polariza- 
tions have been already achieved, and techniques to 
improve upon this have been considered |l0|. Our calcu- 
lations indicate that the dominant source of error may be 
mitigated through narrowing the mixture of dark states. 
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